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Q\ ', A theory determining the electric and magnetic properties of vortex states in Bose- 



> 

o 



Einstein condensates (BECs) is presented. The principal ingredient is the Lagrangian 
of the system which we derive correct to the first order in the atomic centre of mass 
velocity. For the first time using centre of mass coordinates, a gauge transforma- 
tion is performed and relevant relativistic corrections are included. The Lagrangian is 

I symmetric in the electric and magnetic aspects of the problem and includes two key 
Q-i. 

I interaction terms, namely the Aharanov-Casher and the Rontgen interaction terms. 

a , 

^ ' The constitutive relations, which link the electromagnetic fields to the matter fields 
^■ 

via their electric polarisation and magnetisation, follow from the Lagrangian as well 
• I— I . 

^ I as the corresponding Hamiltonian. These relations, together with a generalised Gross- 
ed ' 

Pitaevskii equation, determine the magnetic (electric) monopole charge distributions 
accompanying an order n vortex state when the constituent atoms are characterised by 
an electric dipole (magnetic dipole). Field distributions associated with electric dipole 
active (magnetic dipole active) BECs in a vortex state are evaluated for an infinite- 
and a finite-length cylindrical BEC. The predictd monopole charge distributions, both 
electric and magnetic, automatically satisfy the requirement of global charge neutrality 
and the derivations highlight the exact symmetry between the electric and magnetic 
properties. Order of magnitude estimates of the effects are given for an atomic gas 
BEC, superfiuid helium and a spin-polarised hydrogen BEC. 
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1 Introduction 



Recent years have seen the emergence of the physics of cold atomic ensembles and Bose- 
Einstein condensates (BECs), thanks to advances in laser cooling and trapping [1]. 
Atomic ensembles can now be routinely cooled and confined and, ultimately, made to 
condense under various conditions [1-3], most notably in magnetic and light traps, near 
surfaces, in material cavities and in optical lattices. The trend is towards producing 
larger systems with increasing numbers of constituents and increasing densities than 
currently. Furthermore, continuous evaporative cooling in optically trapped bosonic 
and fermionic gases have been demonstrated in stable optical traps and the techniques 
promise the attainment of BECs with very high densities and much larger coherence 
lifetimes [4]. 

In parallel with these experimental advances, various theoretical investigations [5] 
have been made, with increasing recent emphasis directed towards the angular motion 
of BECs in the form of vortex states [6-9] and the electromagnetic features associated 
with them. In particular, Leonhardt and Piwnicki [7] examined the magnetic effects 
due to the angular motion of a BEC whose constituent atoms are characterised by 
electric dipoles, thus forming a quantum dielectric in angular motion. They predicted 
that such a motion should give rise to a magnetic monopole positioned at the core 
of an n = 1 vortex. More recently, Shevshenko considered the converse case and put 
forward a theory for the generation of vortices in superfluid films by application of 
crossed electric and magnetic fields [9]. 

The situation is now ripe for a rigorous theory of the electromagnetic properties 
of BECs in rotational motion to be constructed. In contemplating such a theory, an 
immediate question which springs to mind is as follows. Hitherto, only the electromag- 
netic properties of dielectric BECs, namely those whose constituent atoms are char- 
acterised by an electric dipole moment (electric BECs), have been considered. What 
if the BEC atoms are characterised by a magnetic dipole moment (magnetic BECs)? 
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This question is now of practical relevance since BECs can be created in light traps, 
allowing, the magnetic properties of, for instance, alkali metal BECs to be investigated 
[8]. More significantly, explicitly magnetic BECs, particularly spin-polarised hydrogen 
BECs, have been created [10]. Intuitively one expects the physics of the problem to 
display a symmetry between the electric and magnetic aspects. Thus if a rotating 
electric BEC generates a magnetic monopole, then a rotating magnetic BEC should 
generate an electric monopole. One of the main tasks of the theory presented here is 
to test the vahdity of such an expectation. 

In this paper we show that rotating magnetic BECs are indeed predicted to give 
rise to electric monopoles and that it is possible to formulate a theory in which both 
types of effect arise in a unified manner. However, our theory predicts that rotating 
BECs should not, as in Ref.[7] give rise to localised monopoles, but, rather, to monopole 
distributions. The treatment automatically ensures that charge neutrality is correctly 
preserved by the electric and magnetic monopole distributions associated with the 
vortex states of the two types of BEC. 

The outline of this paper is as follows. In section 2 we begin by introducing the 
model system as a single atom interacting with the electromagnetic potentials via the 
familiar charge and current sources. We then systematise the procedure based on a 
gauge transformation leading from the conventional Lagrangian of this system to a form 
displaying the convective velocity- dependent interactions. These interactions involve 
the couphng of the electric and magnetic field intensities to the entire series of electric 
and magnetic mutipoles in exact closed forms. In particular, we identify the Rontgen 
interaction [11] and we make use of symmetry arguments to incorporate an Aharanov- 
Casher interaction term [12]. In section 3 we derive the field equations which permit 
the identification of constitutive relations connecting matter fields to electromagnetic 
fields. We also derive the single-atom plus fields Hamiltonian and make use of this 
in section 4 to construct the many-body formulation appropriate for an ensemble of 
such atoms forming a BEC. A main outcome of the many-body theory is a systematic 
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derivation of a generalised Gross-Pitaevskii equation. In section 5 we consider rotating 
BECs and solve the constitutive relations, together with the Gross- Pitaevskii equation, 
to determine distributions for an order n vortex state of the two types of BEG. We 
also display the results for typical situations involving order of magnitude estimates 
of the effects, evaluated for an atomic gas BEG, superfluid helium and spin-polarised 
hydrogen BEG. Section 6 concludes with a summary and provides further comments. 



2 Lagrangian 

It suffices to model the constituent BEG atom in terms of a neutral two-particle sys- 
tem of mass M = mi + m2 involving two bound charges ei = —62 = e of masses 
mi and m2, position vectors Qj and velocities where i — 1,2. The conventional 
Lagrangian of this two-particle system coupled to the electromagnetic scalar potential 
(f) and vector potential A in an arbitrary gauge is written in the following form, which 
is non-relativistic in the motion of both particles [14-17], 

L = Imiql + + / [^60 (E2(r) - c'^^v)) + J(r).A(r) - p(r)0(r)] , (1) 

where E = —A — V0 and B = V x A. The charge and current densities are those 
appropriate for point charges, namely 

= ^i^i^ - q^); = H eiqi(5(r - q^) (2) 

i=l,2 1=1,2 

In order to reveal the bound nature of the two particles and exhibit the desired effects 
involving gross motion of the whole atom we need to express the theory in terms of 
the centre of mass and relative coordinates. We therefore introduce the centre of mass 
coordinate R and internal coordinate q by [16] 

(miqi + m2q2) 

q = qi-q2 (3) 

This step paves the way for the apphcation of a Power- Zienau- Wo oUey (PZW) gauge 
transformation [16] , to be carried out next at the Lagrangian level. As far as the authors 



are aware, this is the first time this method has been performed when using centre of 
mass coordinates, although similar results were obtained using a multi-polar expansion 
on the equivalent Hamiltonian [14]. The generator of the gauge transformation is 
defined by 

X{r) = - fdh'A{r').J^Ur,r') (4) 

e J 

where r') is a gauge vector function given by 

J^mir, r') =e dX(r - R)5(r' - R - A(r - R)) (5) 
Jo 

Note that despite superficial appearances, this vector function is not the electric mul- 
tipolar polarisation field in closed form. The gauge transformation involves replacing 
the vector potential A and scalar potential in the Lagrangian Eq.([^) by A' and 0' 
related to the untransformd potentials by 

A'(r) = A(r) - Vx(r); 0'(r) = 0(r) + ^ (6) 



The transformed Lagrangian L' is 

dxir) 



L' = L- J d^r I J(r).Vx(r) + p(r) 



(7) 



dt J 

It is straightforward, albeit laborious, to show that the terms involving the gradient 
and time derivative of the gauge function are given by 



Vx = A(r) + i 1 dh' eu(r, r') x B(r') 



and 

§ = |A(rO.;^™(r,r')+(R.^)[A(rO.;^™(r,r')] 



with 



-R.A(R) + dV {A(r').:^„(r,r') + R. [(:F^(r,r') - (5k(r,r')) x B(r')]} 

(9) 

eL(r, r') =e dXX{T - R)5(r' - R - A(r - R)) (10) 
Jo 



Substituting for J and p using Eq.(^ and expressing the coordinates qi and q2 in 
terms of R and q using Eq. (^) we have 

L' = -MR^ + -mq^ + f \-eo [e^t) - c^B^r)} + 7^r).E(r) + 

>l(r).B(r) -R.{7?:r) X B(r)}] , (11) 

where m = mim2/M is the reduced mass. The vector fields Vir) and A^r) are, 
respectively, the electric polarisation and magnetisation associated with the charge 
and current sources of the atom, both of which emerge here as closed integral forms 
representing the entire electric and magnetic multipole series to all orders 

nr)=E-^q^>r) (12) 



1=1,2 



^e)(qi,r) + — 0(q2,r) 



X q (13) 



J^v)= -0(q.,r)x(q,-R) 

i=l,2 ^ 

Note that the magnetisation field vector of the two-particle system involves only orbital 
magnetic multipoles, rather than the spin magnetic moment of the particles, but it is 
possible to generalise the theory to incorporate spin. 

The interpretation of each of the various terms in the gauge-transformed La- 
grangian in Eq.(|lT]) is straightforward, except for last term, identified as the Rontgen 



interaction term [7,11], which involves the coupling between the centre of mass motion, 
the magnetic fiux density and the electric polarisation field of the system. 

The presence of the Rontgen interaction term makes the Lagrangian unsymmetric 
between the electric and magnetic properties, specifically as far as the centre of mass 
motion is concerned. This lack of symmetry suggests that an interaction term is miss- 
ing which, like the Rontgen term, should be first order in R and which couples the 
magnetisation to the electric field via the centre of mass motion. The missing term is 
identified as the Aharonov-Casher term [12] and is obtainable by adding a relativistic 
correction to the Lagrangian in Eq.(|l]) so that p p' where 

p'(r) = p(r) + 4R-{V X A^r)}. (14) 
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Including this correction in the original Lagrangian, and following the same procedures 
as descibed above, leads to the appearance of a new term in the form of the Aharonov- 
Casher term so that the new Lagrangian should be 

■1 



L' = ^MR^ + hnq^ + / dh 



eo E^(r)-c^B^(r) +nr).E(r 



+A<r).B(r) - R. T^r) x B(r) - ^A<r) x E(r^ 



(15) 

which is seen to be symmetric between electric and magnetic interactions. Note, how- 
ever, that the Rontgen term was obtained without introduting relativity while the 
Aharonov-Casher term required a relativistic correction. 

3 Single atom plus fields system 
3.1 Canonical momenta 



The canonical variables in the transformed Lagrangian are the centre of mass coordinate 
R and internal coordinate q for the particle system and the vector and scalar potentials 
A and for the fields. The corresponding canonical momenta are as follows. For the 
centre-of-mass motion we have 



df r ( 1 

P = ^ = MR - £r\ 7\r) x B(r) - ^A<r) x E(r) 



while for the internal motion we have 



dL' 
(9q 



mq ■ 



— 0(qi,r) + — 0(q2,r) 



1 . 

Bfr) - — R X E(rl 



For the fields we obtain the canonical momentum correponding to A as 



dC 

aA 



-eoE-V- \r X M 



(16) 



(17) 



where C is the Lagrangian density, identified in L' as the integrand in Eq. (^51) . The 
canonical momentum corresponding to vanishes. 



3.2 Constitutive relations 



The field equations follow from the Euler-Lagrange equations using and A as canon- 
ical variables. They are identified as the familiar Maxwell equations 

OB 

V.D = 0; VxH = — (19) 

dt ^ ' 

where D is the electric displacement field and H is the magnetic field intensity which 

enter Maxwell's equations Eq.(^) provided that the following constitutive relations 

hold 

H(r) = — B(r) - A1+ R X P (20) 

D(r) =eoE(r)+P+^Rx Al (21) 

It is seen that these relations are themselves symmetric between the electric and mag- 
netic contributions and we note, in particular, that the last term in Eq.(^) arises 
directly from the Aharanov-Casher Lagrangian interaction term The corresponding 
term in Eq.(|20|) arises from the Rontgen Lagrangian interaction term. 



3.3 Hamiltonian 



The Hamiltonian of the atomic system interacting with the electromagnetic fields now 
follows from the Lagrangian by use of the canonical prescription 

H = P.R+p.q + n.A-L' (22) 

which yields, after some manipulations, 

where Tiint accounts for the coupling between the electromagnetic fields and the atomic 
system (ignoring diamagnetic terms) 

H,nt = y'rf3r(^^{P(r)}2 + i7^r).n(r)-A<r).B(r)) 

8 



+ ^Jd'r [P. {T^r) X B(r) + ^oA<r) x n(r)} 

+ {T^r) X B(r) + /ioMr) x n(r)} .P] (24) 

where the magnetisation field Al should now be symmetrised in terms of the internal 



coordinate and its canonical momentum. Also in Eq.(^J), pending further clarification, 
the term involving the square of the polarisation is deliberately included in the inter- 
action Hamiltonian, although part of this term accounts for the Coulomb interaction 
between the particles. 



4 Many-body formalism 

4.1 Quantum field theoretic Hamiltonian 

The single atom plus fields Hamiltonian now needs to be generalised to the many-body 
situation involving an ensemble of atoms plus fields interacting quantum-mechanically. 
We begin by introducing the boson field opertor \1/ describing the internal as well as 
the centre-of mass motion of the ensemble, written as a sum over product eigenstates 
as follows 

^(q,R) = V^(R)x(q) = E ^n,(R)Xn,(q)a„„n, (25) 

ni,n2 

where ipni(R) are state functions associated with the centre of mass and Xnzi^) with 
the internal atomic motion. The labels rii and n2 incorporate all quantum numbers 
specifying these states. The operators ani,n2 and aj|,^^„2 annihilation and creation 
opertors satisfying bosonic commutation relations. Completeness demands that we 
must have 

EC,(R)^n.(R') = '5(R-R'); Ex:,(q)xn.(q') =%-q') (26) 

ni n2 

In the next step the many-body Hamiltonian operator is written as follows 
H = I d3q^3j^^t(q,R)[^Zl + ^^^(q,R) + le„|rf3r|-Lif(r) 
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+ nZ, (27) 



+ 2 



where Ti-f^t is given by 

n?„, = I dh 



-^{V{v)Y + -7^r).n(r) - A<r).B(r) 

2eo ^0 



^ y d^^qrf^R y rf3r ^^(q,R) [P. {^^r) x R{v) + iiqMj) x n(r)} 
+ {TKj) X B(r) + /ioA<r) x n(r)} .P] ^(q, R) (2J 



^ 2M 



with the polarisation and magnetisation converted to field theoretic operators as fol- 
lows: 

-p = V / rf3qrf2R^"f(q, R) d\{qi - R)S{r - R - A(qi - R)^(q, R) (29) 
i=i,2 -^0 

M=J2 7r-f rf'qrf'R^^(q, R) [©(q., r) x p - p x 0(q„ r)] ^^(q, R) (30) 

i=l,2 * 

Note that although the polarisation and magnetisation operators are expressed above 
in terms of qi and q2, they can readily be written in terms of the relative and centre-of 
mass coordinates using Eq.(H). 



We have introduced two new terms in the many-body Hamitonian Eg . (|27|) . The 
term involving Vtrap accounts for the trapping potential used to confine the atomic 
ensemble, and the term involving U (R — R') is identified as the hard sphere collision 
term. This is usually taken to be of the form 

U{R-R') = '^^^6iR-R') = UoS(R-R') (31) 

where a is the scattering length. 
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4.2 Generalised Gross-Pitaevskii equation 



Having arrived at the appropriate many-body quantum field theoretical Hamiltonian 
we can now derive the Schrodinger equation satisfied by the atomic field ^(q, R). This 
formally follows from the Heisenberg equation 



H , * 



(32) 



It is instructive to check first what the outcome would be in the absence of all electro- 
magnetic terms plus terms involving the polarisation. In this case Eq.(|32D yields 

p2 p2 



+ ^ + Vtrap + Uo^\q, R)^(q, R) 



^(q,R) = E^(q,R) (33) 



2M 2m 

Apart from the implicit dependence on the internal quantum states x(ct) and the 
explicit appearance of the kinetic energy term p^/2m, Eq.(p3D is reminiscent of the 
Gross-Pitaevskii equation prior to assuming that all ensemble constituents occupy the 
ground quantum state. 

On reintroducing the electromagnetic interaction terms the full Schrodinger equa- 
tion obtained turns out to be a modified Goss-Pitaevskii equation. The modifications 
involve the presence of the internal states, including their coupling to the gross motion 
and the coupling of the two motions to the electromagnetic fields. In order to derive 
this general equation we need to discuss the origins of various terms which will appear 
in it. First there must be a term accounting for the binding of the particles within a 
single atom, i.e. intra-atom Coulomb interactions, and there must also be terms ac- 
counting for the interaction between the atoms, i.e. the inter-atom interactions. Both 
Coulmb interactions arise from the field theoretic term 

f f r^bll 12 I I'h 12 



= — y [[Vl,J' + [V.nterY) (34) 

We have dropped the additional term involving the transverse part of the intra-atom 
polarisation, as this can be shown to lead to infinite self energies. The full inter-atom 
polarisation term gives the electromagnetic interactions between the atoms. 
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Assuming that inter-atomic separations |R — R'| are typically much larger than a 
dipole length |q|, which amounts to the dipole approximation, we can write 

Nd\^l^{R)\^■, M=iVyLi|^(R)|' (35) 

where N is the number of atoms in the ensemble, d is the electric dipole moment and /x 
is the magnetic dipole moment. The modified Gross-Pitaevskii equation in the dipole 
approximation turns out to be in the form 

+ Vtrap + 1 + C/o*^(q, R)*(q, R) + Vm + 



2M 2m 47reog 



*(q,R) = £;*(q,R) 

(36) 

where V^d is an inter- atom interaction term 



Vm = ^ ^'q'(ci.q') {2^t(q'^ K)^(q'^ r) _ ^t(q'^ K)^(q'^ r)} (37) 

while Vgjn includes the coupling of the electric and magnetic dipole systems, pertaining 
to the internal dynamics, with the centre of mass and the electromagnetic fields 



Vem = -d.E^(R)-/u.B(R) 

1 



+ 



2M 



P.{d X B(R) - ^/ix E(R)} -h {d X B(R) - ^^ux E(R)}.P 



(38) 



Next we assume that the number of atoms in the BEG is very large and that they all 
occupy the ground state. Under such circumstances only the ground state operators 
are involved and these become c-numbers such that 

ao,o ~ Vn ; ~ ^ (39) 

All terms involving other operators a„j^„2 ^nin2 where rii 7^ and 712 7^ give 
vanishing results. Although this greatly simplifies the problem, the internal and the 
centre-of-mass motions are still coupled. To achieve a decoupling of the two motions 
we need to perform enemble averaging. This allows us to write 

Vaa = -mum' [ rf'q'd.q' |xo(q')l' (40) 
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The ensemble averaging also results in the substitution 

E^E^^^Sq (41) 

where E'^ is the chemical potential and Eq is the lowest internal energy, i.e. the ground 
state energy eigenvalue of Schrodinger's equation for the internal motion 



Xo(q) = <^oXo(q) (42) 



[ 2m 47reo5 

with Xo(q) the corresponding hydrogenic ground state eigenf unction. The final step is 
to assume that in a condensate there should be a maximum correlation between the 
atoms, in which case we can identify Vaa as an effective dipole-dipole interaction in the 
form 

VM^-N\^Pom\' (43) 
Co 

The modified Gross-Pitaeveskii equation now becomes 

^(q,R) = £;^*(q,R) (44) 



+ Vtrap + N |Vo(R)|' (- + f/ol + 



_2M ""^ '""^ ^' Veo 
Like the Gross-Pitaevskii equation in the absence of electromagnetic effects, this is 
a non-linear Schrodinger equation. The additional terms are V^d, the dipole-dipole 
interaction, and Vem-, the electromagnetic interaction. All the new terms have been 
arrived at systematically within our theory, including the dipole-dipole term which in 
previous considerations has been added phenomenologically. To our knowledge the new 
interaction terms included in Vem have not been derived rigorously before. To see the 
effects of the electromagnetic fields on the properties of a BEG we now consider the 
special case of a BEG in a vortex state 



5 Rotating BECs 

In view of the typical experimental arrangements for generating BEGs we may assume 
that the BEG is confined within a cylinderical region \z\ < zq with cylinder radius Rq 
and height 2zq. The length zq can be larger or smaller than Rq depending on the type of 
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trap, but below we shall consider the cases where zq ^ Rq and infinite. We perform 
a Madelung transformation [6] and, assuming that the wavefunction is independent of 
z, write 



where (r, 0, z) are the cylindrical coordinates of R. On considering solutions of Eq. 
we may drop the electromagnetic field interaction term, which introduces small second 
order corrections [7], and the dipole-dipole interaction term which is also negligible 
compared to the s-wave scattering term. We thus have for the velocity profile of the 
BEG in an order n vortex state 

and also the radial wavefunction satisfies the equation [7] 



Id 2ME , 



m = 0, (47) 

with a the scattering length first encountered in Eq. (|3T|) ; ^ = r/Ro, riid = N/2zq is 
the one-dimensional atomic density parameter and we set iplC, = 1) = 0. We note that 
the solution to this equation for ni^a = would be the Bessel funtion of the first kind, 
Jn{x), with a root at ^ = 1. Figure 1(a) displays the modulas squared of the numerical 
solution of Eq. (|4^) for various values of ni^a. 



5.1 Electric BEC 

In an electric BEC the constituent atoms are characterised by an electric dipole mo- 
ment. We assume that the moment vectors all point in the z-direction. We concentrate 
on the static case and , setting Ai = 0, we take the divergence of both sides in Eq. 
(^). The divergence of the magnetic fiux density B is automatically zero. Hence, 

V.H = V.(R xV) = V.iV X R) - R.(V x V) 



hn 
Mr 
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where we have expressed the right hand side in the first equahty of Eq. ( P8|) as a 



magnetic monopole charge density divided by jj,Q. The magnetic field intensity can be 
written as the gradient of a scalar field: H = — V$, where $ is the magnetic vector 
potential. This is similar to the result obtained in Ref.[7]. However, in [7] it was argued 
that, although the wavefunction vanishes on the cylinder axis, some atoms will leak 
into the core of the vortex and this gives rise to the monopole charge. The variations 
of the wavefunction across the cylinder were ignored. In solving Eq. (|48|) we note that 
the term involving the delta function in the first equality cancels with an identical term 
arising from integrating the derivative which is proportional to [^/^(O)^. This feature of 
our theory guarantees 'charge neutrality', in contrast with the case in Ref. [7]. 

Figure 1(b) displays the variation of the derivative of the modulus squared of the 
wavefunction. This quantity is proportional to the two-dimensional monopole charge 
density in a symmetry plane containing the cylinder axis and can be written in terms 
of ^ = r /_Ro as follows 



Clearly as the scattering length becomes longer or the number of atoms increases the 
wavefunction broadens and the associated charges accumulate near the axis and edge of 
the cylinder. However, the integral of the density over the cylinder volume reduces to 
zero, as demanded by charge neutrality. This magnetic monopole charge distribution 
can be thought of as akin to that found in a cylindrical capacitor with one set of 
charges on an inner thin cylinder and another set of charges of opposite sign on an 
outer cylinder. A useful measure of the effect is the cumulative positive charge, Qm, 



in the inner region of the cylinder, which can be estimated by integrating the charge 
density up to r = Ro/2 (^ = 1/2). We find 




(49) 




(50) 



For an infinite cylinder, Eq. ( ^HD can be solved using Gauss's law to obtain 




(51) 
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where f is a cylindrical radial unit vector. We see that H is proportional to iV'lOP 
which vanishes outside the cylinder. There is a magnetic field only inside the cylinder, 
another signature of the conservation of 'magnetic monopole charge'. 

If we wish to detect the magnetetic effects associated with the vortex, it may be 
necessary to probe the region outside the BEC which can have a magnetic field for a 
cylinder of a finite height 2zo. The solution for the magnetic potential in the finite 
cylinder case can be derived using Green functions [18]. We find 



Jo Jo 



1- z 1+z 

{h'P + (1 - Z)2)V2 + (/l2J2 + (1 + ^)2)1/2 

(52) 



where is a scaling magnetic potential controlling the order of magnitude of the 
effect 

- SnmRl ^^^^ 
the variable z is in units of Zq and we have used the notation 

h'^e+e-m'^os{ct>% (54) 

Figure 2 displays a contour plot of the magnetic potential on a symmetry plane through 
the cylinder axis and above the symmetry plane oi z — As expected, because of 
the finite height of the cylinder, the magnetic potential is not confined to the inside 
but leaks to the outside of the cylinder. Within the cylinder the magnetic potential 
exhibits its largest variations where the wavefunction changes most rapidly. 

It is useful to obtain order of magnitude estimates of the effects just described 
for the cases of a typical atomic gas BEC and for superfiuid hehum. In the case 
of an atomic gas BEC, we choose typical parameters appropriate for ®^Rb [5] in the 
n — 1 vortex state, assuming that the BEC is confined in a cylinder of dimensions 
Rq — zq — 2//m. Each atom is assumed to be characterised by a transition electric 
dipole moment d = eaB, where as = 0.53A is the Bohr radius, and the s-wave 
scattering length is taken to be a = 59A. Setting ni^a = 100 gives a linear density of 
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riid = 1.7 X lO^^m ^ and a total number of atoms N = 7 x lO'^. Thus we have for $^ 

= 3.3 X 10"^^ A (55) 

and for the positive magnetic monopole charge residing in the inner region of the 
cyhnder, Eq. ( pD] ) 

Q^ = 1.3x 10-28|^(^ = l/2)|2Vs (56) 
which is equivalent to a magnetic field of order 10~^^T. 

In the case of superfluid helium, which is taken to be characterised by macroscopic 
electric susceptibility the polarisation field vector can be written as 

no = xnm\' (57) 

where E is an applied electric field and x = 0.052 is the value of the susceptability. 
We have, in effect replaced the volume density times the dipole moment with a sus- 
ceptability times an externally applied electric field. Thus, for every lVm~^ of applied 
field the magnetic monopole charge in the inner region of the cylinder is 

= 1.2 X 10-2^(^ = 1/2) TVs (58) 

and a magnetic field of order lO^^^T (equivalent to the values in [7]). This larger 
magnitude stems from the smaller helium mass and larger density (approximately 
lO^^m"^ in the case of the atomic BEC above and lO^^m"^ for superfluid helium). 
The potential estimated in Eq.(^) for the Rb atomic BEC is clearly too small to 
be measured on the basis of current experimental capabilities. However, one expects 
that denser atomic gas BECs will become available in the future. The corresponding 
estimate found above in the case of superfluid helium indicates that the effects are, 
indeed amenable to experimental detection, as deduced in [7]. 

5.2 Magnetic BEC 

Next we consider to the case of a BEC in which the constituent atoms are charac- 
terised by magnetic dipoles aligned along the 2;-direction. Setting P = 0, we take the 

17 



divergence of both sides in Eq. (0), and following the procedure in the electric BEC 
case, we obtain 

eoV.E = -^V.(Rx M) = fA4(V x R) - R.(V x 



hn 



•-0 



c^Mr 

It can be seen that, for the fields and potentials, Eq.(^) transforms to Eq . (|48|) when 
d is replaced by — yUo/^- Therefore we can deduce the electric field associted with the 
vortex state of an infinitely long cylinderical BEC directly by simple substitution. The 
electrostatic potential is related to the electric field by , E = — V$ and is given directly 
by by Eq.(|5|) and (|3D, with replaced by 

Figure 1(b) now serves to display a function proportional to the areal electric charge 
density in a symmetry plane such that 

Figure 2 should now be taken to present a function proportional to the electrostatic 
potential. The corresponding cumulative negative electric charge residing in the inner 
region of the cylinder obtained by integrating pe over the cylinder volume is approxi- 
mately 

As to the order of magnitude, the best candidate for observing this effect is spin- 
polorised hydrogen [10] which has been produced as a BEC. The hydrogen is doubly 
spin-polarised, thus = 2/ib = he/rrie = 1.9x 10~^^Am^, atypical density is5xl0^^m~^ 
and the s-wave scattering length is a = 0.72A . Again assuming that niaa = 100 and 
Zq = 5mm we obtain 

$^ = -2.7 X 10"^^ V and Qe = -1.9 x lO^^Vl^ = 1/2)|^C (63) 
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Clearly these magnitudes are too small for the effect to be experimentally detectable 
at present. However, the effect should become amenable to measurement whenever 
denser spin-polarised hydrogen BECs can be produced (with densities of the order of 
those of superfluid helium). 
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6 Comments and Conclusions 



We have presented a rigorous treatment of the electromagnetic properties of BECs 
excuting rotational motion when in an order n vortex state. The theory is based on a 
Lagrangian which we have constructed in such a manner as to emphasise the symmetry 
of the problem with regards to the centre of mass motion. The Lagrangian led us to 
constitutive relations which link the fields to the polarisation sources, including the 
convective contributions arising from the centre of mass motion. We have also shown 
how direct generalisation to the many-body situation can be done, leading to a modified 
Gross-Pitaevskii equation which incorporates familiar interactions plus the coupling to 
the electromagnetic fields. Applications of the theory to the case of a rotating BEG 
is straightforward and we have shown how this leads to the electric and magnetic 
properties of a rotating BEG, depending on whether the constituent atoms are electric 
dipole active or magnetic dipole active. 

In particular, we have seen how a dipole-active BEG generates electric and mag- 
netic monopole charge distributions and their associated electric and magnetic fields 
when the condensate is in a vortex state. We have also shown how, even if no conden- 
sate atoms occupy the vortex core, monopole charge distributions, both electric and 
magnetic, arise which are globally neutral. The problem becomes similar to that of 
finding the potential in a cylindrical capacitor. In the case of an atomic gas BEG of 
finite size of either kind, the vortex state generates fields outside it which are small 
when calculated on the basis of current estimates. Their experimental detection will 
have to await the advent of denser atomic gas BEGs than are currently available. The 
predicted magnetic fields are much larger for a denser and more electrically polarisable 
system such as a superfluid. In the case of superfiuid helium, we have shown that the 
magnetic distributions are sufficiently large to be amenable to experimental detection. 
It is envisaged that vortices generated in future denser BEGs of spin-polarised hydrogen 
would give rise to electric distributions which could be measurable. 
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Figure Captions 



Figure 1 

Variations of (a) the modulus squared of the normahsed BEC wavefunction in the n = 1 
vortex state and (b) its derivative with niaa — 100 (sohd curve), rii^a = 10 (dashed 
curve) and riidd = 0.1 (dotted curve). 

Figure 2 

A contour plot within a symmetry plane of the cylinder showing the variations of 
the electric (magnetic) potential of a magnetic (electric) dipole active BEC in the 
n — 1 vortex state which occupies the shaded region. The potential exhibits angular 
symmetry around R = 0, is symmetric about the plane z = and is in units of $o for 
the electric potential and for the magnetic potential (see text). The condensate is 
such that zo/Rq — 1 and niaa — 100. 
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Figure 1 



Boussiakou 




Figure 2 Boussiakou 



